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BEARING-ONLY  TRACKING  ALGORITHMS 


INTRODUCTION 

A  number  of  studies  in  locating  electronic  emitters  have  been  conducted,  and 
algorithms  (given  in  papers  such  as  Refs.  1,  2,  and  3)  have  been  developed  which  provide 
estimates  of  emitter  locations  based  on  bearing  measurements.  For  the  most  part,  these 
algorithms  are  concerned  with  the  problem  of  locating  stationary  emitters,  or  relatively  slow 
moving  emitters,  based  on  measurements  made  by  a  single  moving  platform.  We  will  modify 
two  of  these  algorithms  such  that  a  moving  emitter’s  location  can  be  estimated  by  using 
bearing  measurements  from  two  stationary  sites.  This  was  accomplished  by  thinking  of  the 
stationary  sites  as  a  single  platform  which  was  moving  rapidly  back  and  forth  between  the 
sites.  A  third  algorithm,  the  modified  Kalman  filter  technique,  was  also  developed.  This 
technique  also  uses  bearing  measurements  from  two  sites  to  estimate  the  emitter’s  location 
but  is  not  dependent  on  the  technology  developed  for  locating  emitters  from  a  single 
moving  platform.  All  methods  involve  some  form  of  triangulation. 

The  performances  of  the  three  algorithms  were  compared  by  operating  them  in 
identical  scenarios  using  a  computer  simulation.  Emitter  motion  was  restricted  to  constant- 
velocity  straight-line  motion.  Random  errors  were  introduced  to  the  measurements  by 
selecting  samples  from  a  normal  noise  distribution  derived  from  a  random-number  generator. 
A  standard  deviation  of  1°  was  assumed  for  the  angular-bearing  noise  distribution.  Bearing 
measurements  were  made  every  2  s  from  alternate  sites. 


STATISTICAL  MOVING  AVERAGE  TECHNIQUE 

The  first  technique  to  be  considered  was  developed  from  an  algorithm  presented  in 
Ref.  1.  In  its  original  form  the  algorithm  was  used  by  a  moving  platform  to  estimate  the 
position  of  a  stationary  emitter.  It  was  also  demonstrated  that  the  algorithm  could  be  used 
to  detect  emitter  motion.  However,  no  attempt  was  made  by  the  authors  of  Ref.  1  to  use 
the  algorithm  in  a  tracking  mode,  since  it  provided  erratic  locatfon  estimates  when  the 
speed  of  the  emitter  was  of  the  same  order  as  that  of  the  platform. 

If  one  considers  the  case  of  two  sites  that  are  sharing  bearing  information  on  a  single 
emitter,  one  can  think  of  the  sites  as  a  single  platform  which  is  moving  rapidly  back  and 
forth  between  the  sites.  This  leads  to  the  following  adaptations  of  the  algorithm.  In  Fig.  1 
we  have  two  detecting  sites  at  (*lt  yj )  and  (x2,  >2)-  Without  loss  of  generality  the  sites 
have  been  located  on  the  x  axis.  An  emitter  is  at  (x,  y),  and  the  current  best  estimate  of  its 
position  is  (x,  y ).  The  angle  H  corresponds  to  the  heading  in  the  single  moving  platform 
case  ( 1  ]  and  is  measured  clockwise  from  north.  The  bearing  angle  B,-  is  measured  clockwise 
from  the  baseline,  and  the  angle  (ff  +  B,)  is  measured  clockwise  from  north.  The  true, 
measured,  and  estimated  bearings  are  given  by  B,,  B,,  and  2,-  respectively,  with  the  sub¬ 
script  i  being  incremented  each  time  a  site  makes  a  measurement.  The  error  in  the  bearing 
measurement  is  denoted  by  A Bj,  that  is,  2,-  =  B,-  +  AB,',  and  A B,-  is  assumed  to  be  normally 
distributed  with  zero  mean.  The  difference  between  the  estimated  position  of  the  emitter 
(x,  y)  and  the  true  position  (x,  y)  is  denoted  by  (Ax,  Ay): 
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Fig.  1  —  Geometry  and  nomenclature  for  the  statistical 
moving  average  technique 


and 


A*  *  x  -  x. 
Ay  =  y  -  y. 


To  further  simplify  the  procedure,  the  following  notation  is  adopted: 


rrtj  •  y  -  yjt  m,  =  V  ~  yj, 
n,-  =  X  -  Xj,  nf  »  x  -  Xj, 

where  j  is  the  site  number  of  the  site  making  the  ith  measurement.  From  Fig.  1 


and 


cos  ( H  +  Bj) 


cos  (H  +  B(  ) 

Solving  (5)  and  (6)  for  B,  and  §j  yields 


and 


-H  +  cos’1 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 
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Making  the  substitutions 


Bi  =  8,  -  AB,-, 
x  =  x  -  Ax , 
y  =  y  -  Ay 


in  (7)  gives 


g  -  A B,  =  -H  +  cos’1 


Ay 


v/(m,  -  Ay)2  +  (n,--  Ax)2 


(9) 


(10) 


The  quantity  (g  -  AS,)  can  be  considered  as  a  function  F,(Ax,  Ay)  that  can  be 
expanded  in  a  Taylor’s  series  about  the  point  (0,  0),  giving 


FjiAXj,  Ay,)  =  Fy(0,0)  +  3F,/3Ax| 


Ax=Ay =0 


Ax,  +  bFj/bAy 


|Ax  =  Ay=0 


Ay,-  +  0,-,  (11) 


where  0,-  represents  higher  order  terms  in  the  series.  For  Ax  *  Ay  =  0 

/  *  V 

m,- 


Fj(0,  0)  =  -H  +  cos 


Making  this  substitution  in  (11)  gives 
g  -  AB,  =  Bf  +  3F;/3Ax 


-l  i 


■s/mf  +  n?  / 
Ax,  +  3Fj/3Ay 


=  B,-. 


Ax=Ay=0 


(12) 


Ay  +  0, .  (13) 


|  A* = Ay  -  0 

Taking  the  partials  and  rearranging  (13)  yields 

8  -  8j  =  [-m,(sgn  n,)/(m?  +  nf  )]Ax,-  +  [jn,-l(mf  +  nf)]Ay,-  +  0(  +  AB,. 


(14) 


Assuming  that  0,  «  AB,,  one  can  extend  the  assumption  on  AB,  to  conclude  that 
Ei  =  (0,-  +  AB,-)  is  normally  distributed  with  zero  mean. 

The  substitutions 


Y;  =  g,  -  B,,  xn  =  -rhjHmf  +hf),  and  x,2  =  nflmf+nf)  (15) 


gives  a  set  of  N  equations  for  N  measurements: 

Yx  -  XnAx  +  Xl2Ay  +  Ex  , 
y2  “  x2\Ax  +  ^22 Ay  +  ^2 » 

(16) 

“  -^mAx  ♦  XN2 Ay  +  . 
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In  matrix  notation  Eqs.  (16)  can  be  reduced  to 

Y  =  Xp  +  E ,  (17) 

where 


■*n 

*12‘ 

*r 

y2 

*21 

*22 

'Ax' 

Y  * 

... 

,  x  = 

... 

... 

,  E  = 

... 

,  and  P  = 

Ay. 

.V 

_*N1 

*2V2_ 

_en_ 

In  Eq.  (17)  J3  and  E  are  unknown;  however,  the  assumption  that  E  is  normally  distributed 
allows  an  estimate  of  |3.  The  standard  least-square  estimate  [4]  of  p  is 

3  =  (XrX)_1XTY.  (18) 


Starting  with  an  initial  estimate  (x,  y)  of  the  emitter’s  position  and  a  (Ax,  Ay)  as 
determined  by  Eq.  (18),  one  can  get  a  new  estimate  of  the  emitter’s  position  (5c,  y): 


x  =  x  -  A*, 

This  allows  one  to  recalculate  X  and  Y  from 


A  _  A  A  A 

y  =  y  -  Ay. 


(19) 


Y  =  8,  -  B, ,  rhj  =  y  -  yj 


and 


(20) 


This  process  can  be  repeated  until  some  measure  of  convergence  is  satisfied.  Random  errors 
in  the  measurements  generally  cause  the  process  to  converge  on  a  point  which  is  not  the  true 
emitter  position  but  is  on  the  average  the  best  estimate  achievable  with  the  assumptions  made. 

Up  to  this  point,  emitter  motion  has  not  been  considered;  however,  the  technique  can 
be  expanded  to  consider  moving  emitters  by  hypothesizing  an  “average”  position  of  the 
emitter  during  the  time  interval  in  which  N  measurements  are  taken.  This  is  the  approach 
which  has  been  applied  to  moving  emitters.  In  essence  a  statistical  moving  average  of  the 
emitter’s  position  is  determined.  Initially  a  rough  estimate  is  made  for  the  (N  +  1  )/2  position 
based  on  N  measurements  made  from  two  remote  sites,  and  this  is  used  in  the  algorithm  as 
outlined  for  the  stationary  emitter.  After  the  algorithm  has  converged  to  a  “better” 
estimate,  a  new  measurement  is  taken,  and  the  oldest  of  the  previous  N  measurements  is 
discarded.  The  process  is  then  repeated,  and  the  time  series  of  “better”  estimates  is  used  to 
constitute  a  track. 

A  computer  model  of  the  process  has  been  developed,  and  some  preliminary  results 
are  available.  Figure  2  shows  the  range  and  crossrange  tracking,  errors  for  an  emitter  ap¬ 
proaching  site  2  at  a  speed  of  500  units/s.  The  sites  are  separated  by  20K  units,  and 
measurements  are  made  at  alternate  sites  every  2  s.  The  results  shown  in  Fig.  2  were 
generated  with  N  m  7  and  four  iterations.  Improved  tracking  performance  was  obtained  by 
setting  N  *  15.  This  result  is  shown  in  Fig.  3;  however,  the  price  paid  to  achieve  this  per¬ 
formance  is  an  increase  in  the  time  required  for  track  initiation. 
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Suggested  in  Ref.  1  is  that  the  relative  speed  of  the  target  with  respect  to  the  observer 
might  also  affect  the  performance.  When  two  sites  are  located  as  described  and  alternate 
readings  are  made  every  2  s,  an  effective  observer  speed  of  ±10K  units/s  is  achieved.  To 
investigate  the  effect  of  emitter  speed  on  tracking  performance,  a  scenario  was  designed 
with  an  emitter  moving  parallel  to  the  baseline.  The  spacing  between  the  emitter  path  and 
the  baseline  was  set  at  30K  units. 

Figures  4  and  5  plot  the  errors  in  the  x  and  y  coordinates  of  the  emitter’s  position  for 
emitter  speeds  of  250  and  1000  units/s  respectively.  Contrary  to  the  results  shown  in 
Ref.  1,  tracking  performance  does  not  change  appreciably  for  the  two  speeds  considered. 
This  is  because  with  two  observers  the  triangulation  effect  is  not  significantly  diminished  by 
the  motion  of  the  emitter  if  the  distance  the  emitter  moves  between  measurements  is  small 
with  respect  to  the  length  of  the  baseline.  (For  1000  units/s  and  a  baseline  of  20K  units 
the  emitter  moves  1/10  the  length  of  the  baseline  every  2  s.)  The  track  for  the  fast  emitter 
settles  down  sooner  than  that  of  the  slow  emitter.  This  is  because  the  fast  emitter  moves 
over  a  greater  distance  in  the  measurement  interval;  consequently  some  of  the  measurements 
of  the  fast  emitter  are  made  at  more  accurate  locations  than  all  of  the  measurements  of  the 
slow  emitter. 


MODIFIED  KALMAN  FILTER  TECHNIQUE 

The  second  technique  to  be  considered  is  referred  to  as  the  modified  Kalman  filter 
technique.  Under  this  technique  one  of  the  sites  is  designated  as  the  reporting  site  and  the 
other  as  the  updating  site.  As  shown  in  Fig.  6,  the  measurements  y(k  +  1 )  from  the  report¬ 
ing  site  are  averaged  to  produce  02(ft)-  (Since  the  measurements  are  equally  spaced  in  time, 
02  W  is  coincident  with  01(fe).  If  the  measurements  were  not  equally  spaced,  a  linear 


Fig.  4  —  Positional  error  in  the  x  and  y  coordinates  for  s  crossing  emitter  with  V  •  250  units/s 
when  the  statistical  moving  average  technique  is  used 
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Fig.  5  —  Positional  error  in  the  x  and  y  coordinates  for  a  crossing  emitter  with  V  ~  1000  units/s 
when  the  statistical  moving  average  technique  is  used 


interpolation  procedure  could  be  employed  to  make  02(fe)  coincident  with  /3j(fe).)  j32(fe) 
and  P1(k)  are  then  used  to  triangulate  the  position  of  the  emitter,  and  the  coordinates  of 
the  emitter  are  fed  into  a  Kalman  filter,  which  is  a  tracking  algorithm  that  minimizes  the 
least-square  error  [5] .  The  coordinates  of  the  emitter  corresponding  to  the  ftth  measure¬ 
ment  by  the  updating  site  are  given  by 


?(fe)  =  c(fe)  sin  Pj  (fe)  +  x1  ,  J(k)  =  c(k)cos01(k)  +  yx  ,  (21) 


where 

c(fe)  »  d  cosf32(*)/(sin  (^(fe)  - 

The  state  equation  for  the  Kalman  filter  is 

X(k  +  1)  =  <*>(*)*(*)  +  r(k)A(k),  (22) 
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where 

rjJ 2 
~2 

T 

,  r(fc)  = 

o 
_o 

with  X(k)  being  the  state  vector  at  time  sample  k  consisting  of  position  and  velocity  compo¬ 
nents  x(k),  x{k),  y{k),  and  y{k),  fe  +  1  being  the  next  update  time,  T  being  the  time  between 
updates,  and  ax(k)  and  ay(k)  being  random  accelerations  whose  covariance  matrix  is  Q{k). 
The  observation  equation  is 

Y(fe)  =  M(k)X(k)  +  V(fe) , 

where 

* (fe fi  fi  o  o  oi 

Y(k)  =  ,  M(k )  =  ,  and  Y(fe)  = 

y(fe)J  |_0  0  1  oj  |_ny(fe) 

with  Y(k )  being  the  measurement  at  time  sample  k  consisting  of  measured  positions  x(k) 
and  J(k)  and  with  V(k)  being  zero  mean  noise  whose  covariance  matrix  is  R(k). 

Any  two  successive  measurements  are  not  strictly  independent,  due  to  the  overlap  of 
measurements  from  the  reporting  site  in  the  triangulation  process.  Furthermore  a  bias  can 
be  introduced  by  using  an  arithmetic  average  in  determining  ) ;  however,  these  incon¬ 
sistencies  do  not  seem  to  affect  the  tracking  performance  of  the  algorithm,  and  these 
deficiencies  can  be  removed  if  deemed  necessary. 

The  recursive  algorithm  which  constitutes  the  filter  is  described  in  Ref.  5.  The  elements 
of  the  covariance  matrix 

ax(k)  axy(k) 

*(*)  =  (23) 

0Xy(k)  02(k) 

as  derived  from  Eqs.  (21)  are  (as  shown  in  the  Appendix) 

o\  =  [d/sin2^! -(32)]2  jsin20!  (cos2(31)o|2  +  sin2P2(cos2/32)a|J,  (24) 

a2  =  Jd/sin2  (P1  -  (32)]2 [(cos4 ^1)0^  +  (cos4/?2)a|J  ,  (25) 

and 

axy  =  [d/sin2  d3x  -02)J2{cos3(31(sinp1)o|2  +  cos3(32(sin /32)o|1J,  (26) 

where  the  feth-sample  notation  is  not  shown. 
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The  performance  of  this  algorithm  was  examined  using  the  same  scenarios  used  in  the 
previous  examples.  Figure  7  shows  the  range  and  crossrange  error  for  an  emitter  ap¬ 
proaching  site  2,  and  Figs.  8  and  9  give  the  errors  in  the  emitter’s  x  and  y  coordinates  for 
an  emitter  flying  parallel  to  the  baseline  at  speeds  of  250  and  1000  units/s.  For  both  the 
radial  and  crossing  emitters  the  modified  Kalman  filter  gives  a  significantly  better  track  than 
that  developed  by  the  moving  average  technique. 

EXTENDED  KALMAN  FILTER  TECHNIQUE 

The  third  technique  to  be  considered  is  described  in  Refs.  2  and  3.  In  these  studies 
primary  consideration  was  given  to  the  special  case  in  which  all  the  bearing  measurements 
were  made  by  a  single  moving  platform.  As  described  subsequently,  the  technique  has  been 
applied  to  the  multisite  case  with  two  or  more  sites  making  asynchronous  observations  of  a 
moving  emitter. 


Fig.  8  —  Positional  error  in  the  x  and  y  coordinates  for  a  crossing  emitter  with 
V  =■  250  units/s  when  the  modified  Kalman  filter  is  used 
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Fig.  9  —  Positional  error  in  the  x  and  y  coordinates  for  a  crossing  emitter  with 
V  =  1000  units/s  when  the  modified  Kalman  filter  is  used 


The  geometry  of  the  situation  is  shown  in  Fig.  10.  For  observation  ft  the  bearing 
measurement  (3(fe)  assuming  no  noise  is  given  by 

tan  P(k)  =  [xT(k)  -  x,]/[yr(fe)  -  y,]  ,  (27) 


where  the  target  location  is  given  by  the  coordinates  xT(k)  and  yT(k)  and  the  ith  site  loca¬ 
tion  is  given  by  the  coordinates  jc,  and  y,-.  Equivalently  (27)  can  be  written  as 

[xT(ft)  -  *,]  cos  P(k)  -  [ yr(ft )  -  y,J  sin  0(fc)  =  0.  (28) 

Further  manipulation  yields 

*,cos/J(fe)  -  y,sin/3(ft)  =  xT(k)  cos  fi(k)  -  yT(k)  sin  (3(k) .  (29) 


If  the  left-hand  side  of  (29)  is  considered  as  a  measurement,  it  can  be  rewritten  as 

x(  cos0(fe)  -  y,  sin/3(fe)  =  xr(fe)  cos  P(k)  -  yT(k)  sin  0(fe)  +  N,  (30) 

where  P(k)  is  the  feth  bearing  measurement  corrupted  by  noise  and  N  represents  the  noise  in 
the  measurement  process.  Equation  (30)  can  be  rewritten  in  the  form  of  an  observation 
equation  for  the  Kalman  filter: 

Z(k )  =  M(k)XT(k)  +  N{k),  (31) 

where 


Z(k)  =  Jt,  cos/3(fe)  -  y,  sin/3(fc), 


XT(k) 


xT 

yT 

xT 

yT 
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and 


M(k)  =  [cos  /3(fe),  -sin  0(k),  0,  0]. 

The  corresponding  state  equation,  which  was  described  in  Ref.  5,  together  with  the  observa¬ 
tion  equation  are  in  a  form  resembling  that  of  a  linear  estimation  problem.  Despite  the  non- 
linearities  in  the  observation  equation  an  attempt  was  made  to  apply  linear  estimation 
techniques.  Before  this  was  done,  two  approximations  were  made:  the  matrix  M(k)  was 
replaced  by 

M(k)  =  [cos  0(fc),  -6in  0(k),  0,  0] ,  (32) 

and  the  variance  of  N(k),  denoted  by  ojy(fe),  was  approximated  by  (as  derived  in  the 
Appendix) 

o£(fe)  =  [xfoj  +  o2)  sin2  0  +  (o2.  +  y|jcos2  0 

-2  8mJ(cos0)xiyia^,  (33) 

where  o*(  and  Oy(  are  the  error  variances  in  the  location  of  the  observation  sites  and  o|  is 
the  variance  in  the  bearing  measurement. 

The  measurement  Z(k)  is  essentially  the  perpendicular  distance  from  the  origin  to  the 
bearing  line  (Fig.  10),  and  this  distance  depends  on  the  location  of  the  sites  with  respect  to 
the  origin.  Because  of  this  feature  it  was  advantageous  to  locate  both  receiving  sites  on  one 
of  the  coordinate  axis.  For  sites  on  the  x  axis  yt  =  0,  and  the  variance  (33)  is  reduced  to 

p£(fc)  -  ^c?o|  +  a2^ sin2/)  +  0^  cos2/).  (34) 

With  these  assumptions  the  algorithm  described  in  Ref.  2  was  exercised,  with  results 
as  follows. 
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The  same  scenarios  which  were  described  in  the  two  previous  sections  were  used  to 
examine  the  performance  of  the  algorithm.  Figure  11  shows  the  range  and  crossrange  errors 
for  an  emitter  approaching  site  2  at  a  speed  of  500  units/s.  It  was  assumed  for  all  of  the 
scenarios  that  the  locations  of  both  sites  were  known  almost  exactly;  that  is,  ax.  =  ay.  =  10 
units,  and  a  value  Op  =  1°  was  used  in  computing  the  errors.  Figures  12  and  13  show'the 
errors  in  the  emitter’s  x  and  y  coordinates  for  an  emitter  flying  parallel  to  the  baseline  at 
speed  of  250  and  1000  units/s  respectively.  In  general  the  extended  Kalman  filter  takes 
longer  to  settle  down  than  the  modified  Kalman  filter;  however,  the  track  for  the  slow 
emitter,  moving  parallel  to  the  baseline,  is  more  accurate  once  it  is  fully  developed. 


Fig.  12  —  Positional  error  in  the  x  and  y  coordinates  for  ■  crossing  emitter  with 
V  »  250  units/s  when  the  extended  Kalman  filter  is  used 
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Fig.  13  — Positional  error  in  the  x  and  y  coordinates  for  a  crossing  emitter  with 
V  ■  1000  units/s  when  the  extended  Kalman  filter  is  used 


SUMMARY 

Three  algorithms  have  been  developed  which  use  bearing  measurements  from  two 
remote  sites  to  track  moving  emitters.  Preliminary  results  indicate  that  the  modified 
Kalman  filter  technique  is  superior  to  the  other  two  techniques  in  its  ability  to  rapidly 
develop  a  reasonable  track.  The  extended  Kalman  filter  demonstrated  superior  tracking 
performance  once  it  had  settled  down;  however,  for  the  cases  considered  this  required  an 
unreasonable  amount  of  time.  Further  investigation  is  required  to  assess  the  performance 
of  the  algorithms  against  maneuvering  targets.  The  effect  of  varying  the  time  between 
measurements  should  also  be  considered. 


REFERENCES 

1.  J.  L.  Poirot  and  M.  S.  Smith,  “Moving  Emitter  Classification,’’  IEEE  Transactions  on 
Aerospace  and  Electronic  Systems  AES-12  (No.  2),  255-269  (Mar.  1976). 

2.  A.  G.  Lindgren  and  K.  F.  Gong,  “Position  and  Velocity  Estimation  Via  Bearing  Observa¬ 
tion,”  IEEE  Transactions  on  Aerospace  and  Electronic  Systems  AES-14  (No.  4), 
564-577  (July  1978). 

3.  V.  J.  Aidala,  “Kalman  Filter  Behavior  in  Bearings-Only  Tracking  Applications,”  IEEE 
Transactions  on  Aerospace  and  Electronic  Systems  AES-15  (No.  1),  29-39  (Jan.  1979). 

4.  F.  A.  Graybill,  An  Introduction  to  Linear  Statistical  Models ,  New  York,  McGraw-Hill, 
1961. 

5.  G.  V.  Trunk  and  J.  D.  Wilson,  “Tracking  Filters  for  Multiple-Platform  Radar  Integra¬ 
tion,”  NRL  Report  8087,  Dec.  1976. 


13 


Appendix 

DERIVATIONS  OF  EQUATIONS  (24)  AND  (33) 

DERIVATION  OF  EQ.  (24) 

From  Eq.  (21) 

x(k)  *  csin^  +  Xj 

=  d  cos  02  Pi  /[sin  (/Jj  -  P2)  1  • 

The  differential  form  of  x(k)  is 

Ax  =  d  [~sin/32  sin/3j  /sin  (p1  -  p2)  *  cos/32  sin/3j  cos(0i  ~P2)/sin2  (Pi  -  j32)JA02 
+d  [cos p^  cos  02 /sin  (P1  -p2)  -  cosP2sinp1  cos  (0!  -  02)/sin2  (Pi  ~P2)]AP1  . 
Using  trigonometric  substitutions  and  combining  terms  yields 

Ax  =  d  [sin  j31  cos^/sin2  (Px  -(32)]Aj32  -  d  [sin/?2  cos 02 /sin2  (/31  -02)]  . 

Since  E(Ax)  =  E(A^j)  =  E(  A02)  =  0 

and  a2  =  Z(Ax2),  o^  =  £(40,),  and  Op2  =  E(A02), 

the  result  is 

o2  =  [d/sin2  (0X  -  02)]2  [sin2  P1  (cos2  0X  o^2j+  sin2  02  (cos2  /32jo|J  . 
Equations  (25)  and  (26)  are  derived  similarly. 

DERIVATION  OF  EQ.  (33) 

The  measurement  Zm  can  be  expressed  as 

Zm  -  Z  +  A Z, 

where 

AZ  -  (dZ/dx.JAx,-  +  (3Z/3y,.)Ay,-  +  (3Z/30)A/3 
“  (cos  (3) (A x,-yfA0)  -  (sin  0)(x,A0  +  Ay,) . 
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The  assumption  that  A Z^N  gives 


-2  sin  0  (cos  0)*,yf00  . 
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